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Flaps can be detached from a thin film glued on a solid substrate by tearing and peeling. For flat 
substrates, it has been shown that these flaps spontaneously narrow and collapse in pointy triangular 
shapes. Here we show that various shapes, triangular, elliptic, acuminate or spatulate, can be 
observed for the tears by adjusting the curvature of the substrate. From combined experiments and 
theoretical models, we show that the flap morphology is governed by simple geometric rules. 



Cracks and fractures are commonly observed around 
us in various contexts [1] ranging from drying mud [2] to 
broken windows or ice floes [3J 0] . The classical fracture 
theories, initially formulated by Griffith and Irwin [Sid], 
can reliably predict the onset of cracks motion. In con- 
trast, no general theory is able to predict the path of a 
crack as it propagates. In some specific cases however, 
interesting insight has been gained about cracks trajecto- 
ries. For example, thin films offer a particularly efficient 
set-up to study fracture propagation by limiting the crack 
motion to a two-dimensional manifold. In this context, 
the crucial role of geometry was identified in some oscil- 
latory fracture patterns obtained when a brittle elastic 
thin sheet is cut by a moving object [3 [8]. The propaga- 
tion of two interacting cracks in torn thin films is another 
example where fracture path can be understood [51111]. 

For adhesive thin films, the adhesion strength influ- 
ences the crack paths. It was shown recently how thin 
film elasticity, adhesion and fracture act together to gen- 
erate quasi-perfect triangular tears for pulled adhesive 




FIG. 1: (color online). Typical morphologies of flaps detached 
from adhesive tapes pulled from flat and curved substrates 
with a peeling angle of 90°. (a) Flat substrate, triangular 
shape. When the peeling front propagates by a infinitesimal 
d£ along i-axis, it produces two cracks of length ds together 
with a variation of the width dW and a variation of the peeling 
area dA. (b) Substrate with negative curvature, acuminate, 
(c, d) Substrate with positive curvature with Wo > W c , spat- 
ulate (c) and Wo < W c , elliptic (d), where Wo is the initial 
width of the flap. Negative (positive) curvature, W c , Of and 
<d are defined in the text. Scale bars are 5 mm. The peeling 
front propagates upward. 
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FIG. 2: (color online) . Experimental set-up showing the peel- 
ing/tearing process for positive (a) and negative (b) curva- 
ture. The peeling front advances at a constant velocity v 
with a constant peeling angle a. R is the radius of curvature 
of the peeling front, (c) Cross-section along the cylindrical 
axis of the system, r is the radius of curvature of the ridge 
connecting the flap to the substrate, L is the length of this 
ridge and F is the applied force, (d) Evolution of the tearing 
angle, 8{, for a flat substrate as a function of the peeling an- 
gle, a. Each data point corresponds to the average of at least 
5 experiments. 



tape on flat substrates [T2]. Regardless of the material 
properties, i.e. adhesive/fracture energy, film thickness 
or bending modulus, the triangular morphology is ex- 
tremely robust; only the tearing angle, 9{, is modified. 
This property was recently used to study the elasticity 
of graphene from tearing experiments [13) . 

In this Letter, we show through a combined experimen- 
tal and theoretical study, how geometry of the substrate 
on which the thin film is glued gives exquisite control of 
the shape of the tears. With cylindrical substrates, it is 
possible to generate convergent and even divergent crack 
paths as presented in Fig. [I] 

In the experiments, a thin elastic sheet is adhered 
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to various flat or curved solid substrates. Two paral- 
lel notches are cut on one edge of the adhesive film to 
obtain a rectangular flap. The flap is then pulled at a 
constant speed (v = 1 mm/s) and at a constant peeling 
angle (54° < a < 180°) with a Probe Tack device (see 
Fig. [2^i, b). The final detached flaps were then digitized 
and further analysed (ImageJ and Plot Digitizer) to ob- 
tain the shape of the various tears (selected tears are 
shown in Fig. [lj. 

Flat substrates. The shapes of tears are fully deter- 
mined by the evolution of the flap width, W, with the 
peeled distance, £ during the advancing of both cracks 
(Fig. lj. In order to explain this evolution, we use a 
model based on the Griffith's theory of fracture to de- 
scribe the path of these tears [T2]. It should be noticed 
that due to the peeling geometry, the pulling force in- 
duces a deformation of the adhesive film localized in a 
ridge with a curvature n = 1/r that connects the flap to 
the substrate, see Fig. [2J:. The energy embedded into the 
system by pulling the flap is thus shared among the elas- 
tic energy stored in the ridge, the penalties associated to 
the creation of new surfaces, when cracks advance (7, the 
work of fracture of the film), and the energy dissipated 
in the de-adhesion process (r, the adhesive energy per 
unit area). The variation of energy for an infinitesimal 
motion of the peeling front is, thus, given by 

dU = dU E + 2~ft ds + rdA, 

where Ue, t, ds and dA are the elastic energy of the flap, 
the film thickness, the increment in crack length, ds = 
d£/ cos Of, and the variation of the peeling area dA ^ Wdl 
(see Fig. [J). 

Assuming that Ue depends only on W and £, the vari- 
ation of energy needed to move the front by a distance 
d£ becomes (d x stands for (d/dx)), 
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where dW = -2tan0fd£ since by convention a positive 
tearing angle leads to a decrease of the flap width. To 
solve this equation and find the function 6f(£) that de- 
termines the shape of the flap, we need to compute the 
elastic energy stored in the ridge of curvature k. The av- 
erage curvature of this ridge is simply given by n ~ a/L, 
where L is the length of the ridge (Fig. |2h) . The elastic 
energy is, therefore, given by 



U E ~ BWa 2 /2L 



(2) 



(where B is the bending modulus of the film) . The length 
of the ridge L can be obtained through energy minimiza- 
tion, dU/d£ = 0, reflecting a balance between the different 
energy terms of Eq. (IJ. Considering the expression of 
Ue given above, we obtain the relation 
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Assuming that the width of the flap is much larger than 
the length of the ridge (so that the first term of Eq. ^ 
is negligible compared to the second one) and using the 
approximation cos(9f ~ 1 (valid for the studied range of 
angles), we obtain the expression for L, 
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Inserting this expression in Eq. ^ yields the final ex- 
pression for the elastic energy. 

The propagation of the cracks follows the angle Of 
which minimizes the variation of energy (i.e., this can 
be considered as the maximum-energy-release-rate crite- 
rion |14|). The tearing angle is determined by considering 
that (de f dU) = 0. From Eq. ([I]), this condition becomes 



sin U{ : 



3wU_e_ 
It 



(4) 



The obtained final expression for the elastic energy com- 
bined with Eq. Q determines the tearing angle as a func- 
tion of the material properties 
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Depending on the relative importance of adhesion and 
crack energy terms (determined by the ratio rW/ft), 
we can distinguish two different regimes. For large 
flaps, adhesion dominates the evolution of the crack path 
(W » yt/r - 1 mm in our experiments). This regime 
characterized by a constant tearing angle Of, is described 
by the relation 



sin Of ■ 



2yt 



(6) 



This relation is similar to the one obtained in Ref. [12j . 
the influence of the peeling angle a being added. The 
linear relation between sin#f and a is in close agreement 
with the experimental data shown in Fig. |2h. 

In addition, a regime dominated by the fracture en- 
ergy should always be observed near the tip of the flaps, 
i.e., when W becomes very small. This regime is char- 
acterized by sin^f = a^/B/ytW. The tearing angle thus 
increases for decreasing flap width; this yields a short- 
ening of the flap and a re-entrant morphology (see for 
instance Fig. 1 of [12] ). Considering that sin6*f ^ tan#f = 
d(W/2)/d£, the shape of the tears in this regime corre- 
sponds to the power law W ~ £ 2 ^ 3 . However, it should be 
mentioned that this regime is rather difficult to observe 
when conventional adhesive tapes are used. 

Curved substrates. When a thin film is glued on a 
curved substrate, the peeling/tearing process produces 
various morphologies for the flaps as depicted in Fig. [T] 
Those shapes depend on the radius of curvature, R, of 
the substrate. For clarity, we define the curvature as 



shown in Fig. [2^,b and use, in the following, the param- 
eter e equals +1 or -1 for positive or negative curvature, 
respectively. The main ingredient determining the global 
shape of the tears is the morphology of the peeling front 
induced by the substrate curvature. This front is linear 
for flat substrates, but, in contrast, it is characterized by 
a finite radius of curvature, R, for cylindrical substrates 
as shown in Fig. [3]:. The direction of the curvature with 
respect to direction of the cracks propagation is deter- 
mined by the sign of the substrate curvature. This is 
the main source for the breaking of symmetry which ex- 
plains the occurrence of closing or opening cracks path 
when the peeling of the adhesive film is performed on 
the exterior (e = -1) or the interior (e = +1) side of the 
cylinder, respectively. 

The relation between the radius of curvature of the 
peeling front and the radius of curvature of the substrate 
can be explained with the help of simple geometrical ar- 
guments. The limiting zone of contact between the adhe- 
sive film peeled with an angle a and the substrate is given 
by the intersection between the cylindrical substrate and 
a plane forming an angle a/2 with the symmetry axis 
of the cylinder as shown on Fig. ^p. The projection of 
this curve on the (x, y) plane is characterized by a radius 
of curvature at its tip, R, given by R = R tan a/2, see 
Fig. [3];. This relation corresponds to the geodesic cur- 
vature of the intersection curve |15j and is in very good 
agreement with the evolution of R with the peeling angle 
(Fig. [3jl) for various values of R. The curvature of the 
peeling front is not significantly influenced by the speed 
at which the front advances, i.e. by the adhesion energy, 
indicating that the origin of this curved morphology is 
purely geometric. It is worth noting that the peeling 
front becomes linear for peeling angle a close to 180° ir- 
respective of the curvature of the substrate (R tends to 
infinity for a close to tt). 

In order to derive a model describing the shape of the 
flaps for curved substrates, we make the following as- 
sumptions: i) the tears are produced in the adhesion 
regime. This approximation is justified since almost ev- 
erywhere the width of the flaps is large enough to neglect 
fracture energy (W » 1 mm) except for the extreme tip 
of the convergent tears; ii) the distortion of the peeling 
front induces a rotation of the path direction of the cracks 
by an angle (3 (see Fig. 3^) , keeping the same value of the 
tearing angle 9f, given by Eq. To account for this ro- 
tation of the crack path with respect to the cylinder axis, 
we introduce a new angle O = 0f~e/3 (with e = ±1, see here 
above). The angle /3 is determined by the distortion of 
the peeling front and the width of the flap and is given by 
W/2R. This new geometry transforms the equation that 
determines the shape of the tears for curved substrates 
into 
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FIG. 3: (color online), (a) Origami model of the morphology 
of a thin film adhering onto a cylinder with a peeling angle 
a. (b) Lateral view showing the plane containing the limiting 
zone of contact between the adhesive film and the cylindrical 
substrate and forming an angle a/2 with rc-axis. (c) Top 
view showing the curved morphology of the peeling front. The 
angle /3 defines the rotation at the cracks of the tangent of the 
curved front (blue solid line) with respect to a planar peeling 
front (red solid line). Scale bar: 3 cm. (d) Evolution of the 
observed radius of curvature R/R with the peeling angle a. 



I = £/R, we get the following ordinary differential equa- 
tion (ODE) 
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The solution of this ODE is given by 

W = e[l + (eW - l)e eM ] 



(7) 



(8) 



with W/R as the lowest order and using as above cos#f ~ 
1. Introducing the normalization, W = W/2Rsm9{ and 



where Wo = W(£o) is an arbitrary point along the flap 
profile. 

Let us now discuss in detail the morphology of these 
solutions for the different values of e and Wo- For sub- 
strates with negative curvature, we always observe clos- 
ing tears with an "acuminate leaf-like" shape [T5] and 
a sharp decrease of O with I (Fig. [TJd). With e = -1 
and VF(0) = (i.e., this point represents the ultimate 
tip of the tear), Eq. (8 1 simplifies to W = ~ l] which 
captures almost perfectly the evolution of the tears shape 
for various samples (see the master curve in Fig. [4^,) . For 
the positive curvature, pulling an adhesive can produce 
either opening ( "spatulate" ) or closing ("elliptic") tears 
as shown in Fig. [Tj: and d, depending on the initial value 
of the flap width. Interestingly, these morphologies and 
the transition in shape with the width are also captured 
by the model. For e = +1 , the solution reads 
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FIG. 4: (color online). Master curves of the tear morphologies 
for (a) negative curvature obtained from 45 experiments (R = 
54.8 , 20.2 and 16 mm) and (b) positive curvature obtained 
from 30 experiments (R = 18.5 and 14.1 mm). The scalings 
used to normalize the plots are discussed in the text. The 
solid lines correspond to Eq. Q obtained from the model 
with y = W/2. 

From this equation, two different regimes are indeed ex- 
pected depending on the flap width. For large initial 
width (W > 1 or W > W c = 2_Rsin# f ), we observe a 
rapid opening of the tears contrasting with the trian- 
gular shapes previously reported. This behavior is pre- 
dicted by the model since in this case, the exponential 
term of Eq. Q is strictly positive. The agreement with 
the experimental data is very good, a master curve of 



all the results being given in Fig. [4Jd. For small width 
(W < 1 or W < W c ), we recover closing tears but with 
an "elliptic leaf-like" shape instead of an acuminate one. 
This behavior is also predicted by our model since in this 
case, the exponential term of Eq. Q is strictly negative. 
Again, the model based on geometry arguments captures 
the variety of the tear morphology observed for curved 
substrates. The transition between opening and closing 
crack paths predicted by the model is also experimen- 
tally observed and perfectly reproduced by the theory 
(Fig. [4)d). It should be emphasized that, for positively 
curved substrates and a given set of experimental con- 
ditions (R, a, 7, r, and t) there is a single value of the 
width W c = 2R sin 9{ that produces tears with a constant 
width and parallel cracks. This value obviously corre- 
sponds to an unstable state lying between both opening 
and closing domains. 

As shown here above, the peeling front is linear for 
peeling angle close to tt regardless of the curvature of the 
substrate. As expected from the proposed model, the 
tears adopt a triangular morphology for various curved 
substrates when the flap is pulled at 180°, see the inset 
of Fig. gji. 

In conclusion, we have demonstrated for flaps detached 
from films glued to simple circular cylinders that a con- 
stant mean curvature strongly affects the crack path and 
thus the shape of the tears. Moreover, the proposed 
model can be readily extended to any developable sur- 
face, i.e. conical or cylindrical. Indeed, Eq. Q de- 
scribes the local behavior of the crack propagation and, 
therefore, a spatial variation of the surface curvature can 
be easily introduced. A possible extension of this work 
could be the study of the influence of Gaussian curvature 
on crack path and tear morphology for surfaces such as 
sphere, hyperboloid. 
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Supplemental Material for "How Geometry Con- 
trols the Tearing of Adhesive Thin Films on 
Curved Surfaces" 



I. PEELING FRONT CURVATURE 

In Fig. [5^, below (see also Fig. 3 of the main text), a 
sheet of paper is wrapped around a cylinder to model the 
shape adopted by the thin adhesive film during our peel- 
ing/tearing experiments. The blue solid line indicates 
the peeling front whose shape is identified to a curve 
7 obtained by the intersection between a cylinder and a 
plane forming an angle a/2 with the x-axis (see Fig.[5]D,c). 
When the sheet is unrolled and flattened, the creasing of 
the paper sheet reveals the curvature of the peeling front 
relevant for our study (see Fig. [5ji,e). This procedure al- 
lows to distinguish the part of the curvature of the curve 
7 that is due to the cylinder and the one that comes 
from the peeling front itself. In mathematical terms, this 
means that we are dealing with the geodesic curvature, 
Kg = (Rg)^ 1 , of the curve 7 and not with the total cur- 
vature. As illustrated in Fig. [5^, this geodesic curvature 
varies slowly along the peeling front and is approximated 
in our study by the curvature, i? -1 , at the tip of the front 
(at y = 0). The deviation between the peeling front and 
an arc of circle is also shown. The accuracy of this ap- 
proximation improves as the ratio between the width of 
the front, W, and the radius of the cylinder, R, decreases. 



A. Equation for the intersection curve 7 



The equation of the cylinder is given by 
z 2 + y 2 = R 2 , 



(10) 



with -00 < x < 00. The equation of the plane reads 

z = tan(a/2):r, (11) 

with -00 < y < 00. Consequently, the parametric equa- 
tions for the intersecting curve 7 is found to be 

7(i) = (t,±^R 2 -tan 2 (a/2)i 2 ,tan(a/2)tj , (12) 

with -R/ tan(a/2) <t<R/ tan(a/2). 

B. Expression of R: simple derivation 

Since we are only interested by the geodesic curvature 
at the tip of the curve 7, there are two equivalent ways 
to compute this quantity. The first one is obviously to 
directly compute the geodesic curvature along the curve 
7 by using the general definition of this quantity, as de- 
tailed in the next section. The second equivalent one is 
to compute the curvature at the tip of the curve j xy ob- 
tained by projecting the curve 7 in the plane (x,y), as 
mentioned in the main text. 



From Eqs. (10) and (11), this curve, j xy , is an ellipse 
whose equation reads 



y 2 + tan 2 (a/2)x 2 = R 2 . 



(13) 



The radius of curvature, R, at the tip y = is thus given 
by 



R = i?tan(a/2). 



(14) 



This corresponds to the expression used in the main text 
and tested experimentally (see Fig. 3d of the main text). 



C. Expression of R: geodesic curvature 



An alternative procedure to obtain the result (14) is 



to compute the geodesic curvature defined by the rela- 
tion jT7j 



Kg(t) 
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where prime denotes the derivative with respect to the 
parameter t, \\x\\ is the norm of x and n is the unit normal 
to the cylinder. The radius of curvature R is given by 



R=(\ Kg \r 



(16) 



where n g is evaluated at t = R/ tan(a/2) which corre- 
spond to y = from Eq. ( 12 ). 



The unit normal to the cylinder is computed from the 
gradient of G(x, y, z) = z 2 + y 2 - R 2 as follow 



n(x,y,z) 
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The expression of this normal vector along the curve 7 is 
obtained by replacing x, y and z in terms of t using the 
parametrization ( 12 ) : 



n(t) = (o, a/1 - tan 2 (a/2)(r/ J R) 2 , (t/R) tan(a/2) j . (18) 

The rest of computation is done without difficulty and 
we obtain 

R 9 (t) = K(i)) -1 (19) 
> (l + tan 2 (a/2) -(t/i?) 2 tan 2 (a/2)) 3/2 



R- 



(t/R) tan 3 (a/2) 



Finally, the geodesic curvature at the tip of the curve 7 
is given by 



R = R g (R/ tan(a/2)) = i?tan(a/2) 



(20) 
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FIG. 5: (a) Origami model of the morphology of a thin film adhering onto a cylinder with a peeling angle a. Scale bar: 3 
cm. (b) Lateral view showing the plane containing the limiting zone of contact between the adhesive film and the cylindrical 
substrate and forming an angle a/2 with 2-axis. (c) Intersection curve 7 between a cylinder and a tilted plane, (d) View of 
the wrapped paper sheet around a cylinder (total curvature), (e) View of the flattened paper sheet (geodesic curvature), the 
creasing of the paper sheet reveals the intersection curve which is compared to an arc of circle (blue dotted line). 



